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^3 We propose an extraction of the running coupling constant of QCD in the infrared 

Oh region from experimental data on deep inelastic inclusive scattering at Bjorken x — > 1. 

We first attempt a perturbative fit of the data that extends NLO PQCD evolution to 
large x values and final state invariant mass, W, in the resonance region. We include 
both target mass corrections and large x resummation effects. These effects are of order 
0(1/Q 2 ), and they improve the agreement with the Q 2 dependence of the data. Standard 
analyses require the presence of additional power corrections, or dynamical higher twists, 
— i to achieve a fully quantitative fit. Our analysis, however, is regulated by the value of the 

strong coupling in the infrared region that enters through large x resummation effects, 
V^*) and that can suppress, or absorb, higher twist effects. Large x data therefore indirectly 

V^*) provide a measurement of this quantity that can be compared to extractions from other 

£S^J observablcs. 
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1. Introduction 



The main goal of QCD is to describe the structure of hadrons in terms of its funda- 
mental degrees of freedom given by the quarks and gluons, or partons. Hadrons are 
observable in both the initial and final stages of a hard process while the existence 
and properties of partons are inferred only indirectly. Because of the smallness of 
the running coupling constant of QCD at large enough momentum transfer squared, 
Q 2 , or equivalently at short distances 0(1/ y/Q 2 ), a hard probe sees hadrons as com- 
posed of "free" quarks and gluons carrying fractions x of the hadron's Light Cone 
(LC) momentum, with given probability distributions. Factorization theorems reg- 
ulate this fundamental property of the theory by allowing for the short distance 
behavior at a given scale, Q 2 , to be evaluated separately from the long distance 
contribution which is identified with the Parton Distribution Functions (PDFs) . ^. 
As increasingly shorter distances are probed the PDFs shape in x changes due to 
radiative processes, according to a pattern which is calculable with high accuracy 
within Perturbative QCD (PQCD) E. 

Although the perturbative phases of a hard collision are clearly distinct from 
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the non-perturbative stage that gives rise to hadron structure, experimental obser- 
vations indicate that such non-perturbative stage appears to be almost ubiquitously 
pre-determined by - or to be keeping track of - the perturbative mechanism. This 
concept goes under the name of parton-hadron duality. Parton-hadron duality im- 
plies in the most extreme case, that hadronic observables are replaced by calculable 
partonic ones with little more going into the hadronic formation phase of the pro- 
cess. In this respect, duality can be seen as a phenomenological manifestation of 
the non-perturbative to perturbative transition in QCD. A full explanation for its 
onset is however still lacking. 

As an example in this contribution we will focus on Bloom Gilman (BG) dual- 
ity^, observed in DIS. Here the relevant kinematical variables are: x = Q 2 /2Mv 
(M being the proton mass and v the energy transfer in the lab system), the four- 
momentum transfer, Q 2 , and the invariant mass for the proton, P, and virtual 
photon, q, system, W 2 = (P + q) 2 (W 2 = Q 2 (l/x — 1) + M 2 ). For large values 
of Bjorken x > 0.5, and Q 2 in the multi-GeV region, one has W 2 < 5 GeV 2 , i.e. 
the cross section is dominated by resonance formation. While it is impossible to 
reconstruct the detailed structure of the proton's resonances, these remarkably fol- 
low the PQCD predictions when averaged over x (see e.g. ^for a review). Although 
BG duality was observed at the inception of QCD, quantitative analyses could be 
attempted only more recently, having at disposal the extensive, high precision data 
from Jefferson Lab^. Such analyses fall into two main categories: 

i) Quark model/symmetries-based, RefsHUZIH These analyze the leading order ma- 
trix elements for the quark proton scattering subprocess and study the conditions 
under which the DIS cross section which is proportional to the incoherent sum of 
the squares of the constituent charges, can in average connect smoothly to the res- 
onance production portion of the cross section which is given by the square of the 
coherent sum of the constituent charges. This connection is regulated by Parity 
transformations. It was, in fact, shown how the excitation of resonance states of 
opposite Parity interferes destructively in all cases except for the leading twist. No 
account of perturbative processes enters this model. 

ii) PQCD-based, Refs. I 9 l 10 l n [ This approach claims that only after a complete 
perturbative QCD analysis is performed can one define duality by quantitatively 
establishing whether, and to what extent, this phenomenon is responsible for the 
apparent cancellation of multiparton correlations. It was noticed that the extension 
of NLO PQCD evolution to large x must include large x resummation effects. The 
overall effect of these is to shift the scale at which as is calculated to lower values, 
with increasing x. Standard analyses require the presence of additional power cor- 
rections, or dynamical higher twists, to achieve a fully quantitative fit. Our analysis, 
however, being regulated by the value of the QCD running coupling in the infrared 
region that enters through large x resummation effects, brings to a suppression of 
higher twist effects. These get, in fact, absorbed in the coupling's infrared behavior. 
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Although these two approaches have been considered so far complementary to 
each other, a unified description might derive through the definition of the effective 
coupling 

This paper is organized as follows. In Section 2 we present an overview of the 
large x data, and discuss a few aspects of the evolution mechanism for DIS at large x 
where two scales related to the invariant mass and to the four-momentum transfer, 
are simultaneously present. We then illustrate in Section 3 the connection between 
large x data and the coupling in the infrared region: we argue that once the range of 
validity of parton-hadron duality is defined quantitatively, a precise PQCD analysis 
at large x would open up the possibility of extracting the strong coupling constant 
at low scale. In Section 4 we draw our conclusions 



2. Analysis of large x data 

Our studies of the physical origin of duality, and of its impact on our understanding 
of the nucleon's structure started a few years ago when we set up a program to 
quantitatively extract the scale dependence of the large x Jlab Hall C data^. In the 
analysis performed in Refs 121121 we addressed several effects that have a large impact 
at large x, namely Target Mass Corrections (TMCs), large x resummation effects, 
and higher twists. This work was then completed, and extended to polarized data in 
Ref. An important point emerged from Refs. IHLL2LLU that a deeper understanding 
was needed of those aspects unique to the large x perturbative QCD analysis. 

High precision inclusive unpolarized electron-nucleon scattering data on both 
hydrogen and deuterium targets from Jefferson Lab are available to date in the 
large x, multi-GeV regime (see ^ and references therein). Because of the precision 
of the data one should now be able to distinguish among different sources of scaling 
violations affecting the structure functions in addition to standard NLO evolution, 



• Target Mass Corrections (TMC), 

• Large x Resummation Effects (LxR) 

• Nuclear Effects 

• Dynamical Higher Twists (HTs), 

• Impact of Next-to-Next-to-Leading-Order (NNLO) perturbative evolution. 



All of the effects above can be extracted with an associated theoretical error. It is in 
fact well known that their evaluation is model dependent. Recent studies, however, 
have been directed at determining more precisely both the origin and size of the 
associated theoretical error. Recent analyses have been taking into account, so far, 
some but not all of the effects listed above 
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2.1. Unpolarized structure function 



The inclusive DIS cross section of unpolarized electrons off an unpolarized proton 
is written in terms of the two structure functions F-i and F%, 



(Mxy) 



i y - " Q2 ) l : ~2 + !f-rl- 



(1) 



dxdy Q 2 xy 

with y — f/ex, e\ being the initial electron energy. The structure functions are 
related by the equation 

F 1 =F 2 {l + 1 2 )/{2x{l + R)), (2) 

where 7 2 = 4M 2 x 2 /Q 2 , and R is ratio of the longitudinal to transverse virtual 
photo-absorption cross sections. In QCD, F 2 is expanded in series of inverse powers 
of Q 2 , obtained by ordering the matrix elements in the DIS process by increasing 
twist r, which is equal to their dimension minus spin 

F 2 (x, Q 2 ) = Ff T (x, Q 2 ) + ^ + (i,) * FF(P, Q 2 ) ( 1 + q ^) +0 (^ 

The first term is the leading twist (LT), with r = 2. The terms of order 1/Q T ~ 2 , 
t > 4, in Eq. ^ are the higher order terms, generally referred to as higher twists ^. 

Target Mass Corrections 

TMCs are included in F 2 LT . For Q 2 > 1 GeV 2 , TMCs can be taken into account 
through the following expansion^ 

F^ TMC \x,Q 2 ) = ^F™^Q 2 ) + ^ f (4) 

where F£° is the structure function in the absence of TMCs. A more recent analysis 
023 re-examined TMCs within the collinear factorization approach of^in order to 
address the longstanding question of the unphysical behavior in the threshold region 
of Eq.Q. This originates from the fact that as x — > 1, one obtains a Q 2 -dependent 
threshold, namely F 2 (^,Q 2 ) = for £ > £ max = 2/1 + yj\ +4M 2 /Q 2 , therefore 
rendering F 2 undefinable as Q 2 varies (see discussion in 1211). In the formalism of 
Ref. TMCs, applied to the helicity dependent structure functions read 



F^™ c \x,Q 2 )= * dm 2 P (m 2 )F? 



(5) 



where Ft = F\. The final quark is assumed to hadronize into a jet of mass m.j with 
a a process dependent distribution/smearing function p(m 2 ). In our extraction we 
take the perspective that the evaluation of TMCs is always associated with the 
evaluation of HTs - TMCs should in principle be applied also to HTs - in an 
inseparable way. Therefore we consistently keep terms of 0(1/Q 4 ) I 11 ! 22 !, whether in 
the formalism/prescription of Ref. ^ or of Ref. H(x, Q 2 ), then, represents the 
"genuine" HT correction that involves interactions between the struck parton and 
the spectators or, formally, multi-parton correlation functions. 
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Threshold Resummation 

In order to understand the nature of the remaining Q 2 dependence that cannot 
be described by NLO pQCD evolution, we also include the effect of threshold re- 
summation, or Large x Resummation (LxR). LxR effects arise formally from terms 
containing powers of ln(l — z), z being the longitudinal variable in the evolution 
equations, that are present in the Wilson coefficient functions C(z). Below we write 
schematically how the latter relate the parton distributions to e.g. the structure 
function F%, 

F 2 LT (x, Q 2 ) = |i I dz C W l( x / z > Q 2 )' (6) 

where we have considered only the non-singlet (NS) contribution to F% since only 
valence quarks distributions are relevant in our kinematics. The logarithmic terms 
in C(z) become very large at large x, and they need to be resummed to all orders 
in as- Resummation was first introduced by linking this issue to the definition of 
the correct kinematical variable that determines the phase space for the radiation 
of gluons at large x. This was found to be W 2 = Q 2 (l - z)/z, instead of Q 2 E3H 
As a result, the argument of the strong coupling constant becomes z-dependent: 
a s(Q 2 ) — > Ois(Q 2 {l — z )/ z ) l 25 l 26 l In this procedure, however, an ambiguity is in- 
troduced, related to the need of continuing the value of as for low values of its 
argument, i.e. for z — > 1 * 27 . Although on one side, the size of this ambiguity could 
be of the same order of the HT corrections and, therefore, a source of theoretical 
error, on the other by performing an accurate analysis such as the one proposed 
here, one can extract as for values of the scale in the infrared region. We address 
this point in more detail in the next Section. 

Nuclear Effects 

Theoretical uncertainties in the deuteron are taken routinely into account, and 
are expected to be in sufficient control (see ^ and references therein). Uncertainties 
arise mainly from 

i) Different models of the so called nuclear EMC effect; 

ii) Different deuteron wave functions derived from currently available NN potentials, 
giving rise to different amounts of high momentum components; 

Hi) The interplay between nuclcon off-shellness and TMC in nuclei. 

Finally we did not consider NNLO calculations, these are not expected to alter 
substantially our extraction since, differently from what seen originally in the case 
of F3, these have been proven to give a relatively small contribution to F 2 . 

Once all of the above effects have been subtracted from the data, and assuming 
the validity of the twist expansion, Eq.Q in this region, one can interpret more re- 
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Fig. 1. HT coefficients extracted in the resonance region according to the procedure described 
in the text. HT extracted with only the NLO calculation (squares); the effect subtracting TMC 
(open circles); the effect of subtracting both TMC and LxR (triangles). Shown for comparison are 
the values obtained from the coefficient H obtained in Ref. B2J using DIS data and including the 
effect of TMC. Adapted from Ref. EH 



liably any remaining discrepancy in terms of HTs. Since we extend our x-dependent 
analysis to the resonance region we consider the following integrated quantities 

I rcs ((x),Q 2 ) = T maX Fr(x, Q 2 ) dx (7) 

** ^min 

where is evaluated using the experimental data in the resonance region. For 
each Q 2 value: x min = Q 2 /(Q 2 + W 2 ^ - M 2 ), and x max = Q 2 /(Q 2 + W 2 nin ~ M 2 ), 
where W m - ln and W max delimit the resonance region, and (x) is the average value of 
x for each kinematics. This procedure replaces a strict point by point in x, analysis. 

Typical results from the analysis outlined above are plotted in Figure 1 where 
we show the HT coefficient defined from Eq.([3]) as 

R LT = C(x) = Q 2 [F 2 (x, Q 2 )/F^ T (x, Q 2 ) - l] 

The error in the figure is from the experimental data. No theoretical uncertainty was 
included. However, our results clearly show that the combined effects of TMCs and 
LxR substantially reduce C{x). We take this as illustrative of the accomplishments 
one can expect from the analysis we suggest in this contribution. Essential features 
that emerge are the interplay between the values of «5(M|) and the HTs, the 
relevance of TMCs, and, most importantly, the need to define as in the infrared 
region. All of these features can affect the central values of the HTs reported in Fig. 
1. 
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Fig. 2. Phase space in the evolution of the NS component. The dotted lines are for fey MAX = Q 2 , 
for Q 2 = 10, 100 GeV 2 . The full lines represent the upper limit fcj, MAX = Q (1 — z)/z for 
Q 2 = 10, 100 GeV 2 . In this case, one can see a clear reduction of the allowed fey at large z. 



3. Strong coupling constant at high x 

We now discuss in more detail the working of threshold resummation, and its pos- 
sible impact on the analysis of F% at large x^. Starting from NLO, the coefficient, 
C(z) in Eq.|6| is dominated at large x by terms proportional to [as{Q 2 ) hi(l — z)] n 
which need to be resummed in the perturbative series. The physical origin of these 
terms is in the phase space for the contribution of gluons emission to evolution, 
which become soft as x — > 1. A mismatch in the cancellation with the virtual glu- 
ons contributions ensues. If, however, one carefully evaluates the kinematics for 
gluon emission at large x within a quark-parton model view, one obtains ^3 

q (x,Q 2 )= q (x,Q 2 ) + jT ^ aS ^ 1 £~P qq ( X z ,a s )q(z,kl), (8) 

where Q 2 is an arbitrary initial scale, and W 2 — Q 2 (l — z)/z is the maximum k\ 
in the virtual photon-quark center of mass system, appearing in the ladder graphs 
that define the leading log result. The resulting phase space is shown for different 



Q 2 values in Fig. 2.1 The reduction of the allowed k\ results in a simultaneous shift 
in the argument of as — > &s{Q 2 / (1 — z )/ z ), an d a cancellation of the as(Q 2 ) hi(l — 
z) divergence in the NLO coefficient function. As a consequence of rescaling the 
argument of as one has to consider its continuation into the infrared region l 26 l n l 
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Fig. 3. Left panel: as/ir extracted from the analysis of the large x data discussed in the text, and 
plotted vs. z, Eq.|6j| (upper panel), and vs. W = y/Q 2 (l — z)/z (lower panel). For comparison we 
show the extraction from Ref.ES using Jefferson Lab data at Q 2 =0.7 — 1.1 GeV 2 . 



The left panels of Fig. |2.1| display the results for as used in our analysis for different 
values of Q 2 . We also show, on the right, the extracted value of the effective as from 
the GDH sum rule. We therefore suggest large x evolution in DIS as yet another 
way of defining an effective coupling constant at low values of the scale. A more 
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quantitative analysis to relate different types of measurements, and to study in 
depth the possible process dependence of as is in progress^. 

4. Conclusions 

In conclusion, we believe there is a much richer structure to the scale dependence 
of the nucleon's distribution functions that persists behind the apparent cancella- 
tion among higher twist terms. We started uncovering this structure in the initial 
work of RefsEEH Our analysis opens up the possibility of extracting the effective 
strong coupling at low scale, and to connect it to extractions from different pro- 
cesses, namely Refsl 14 ! 15 ! and RefP. While on one side this points at the fact that 
PQCD provides an essential framework for understanding the working of duality, a 
possible connection with the model of Refs03HHis envisaged through the definition 
of effective charge 
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